Pair production with neutrinos in an intense background magnetic field 
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We present a detailed calculation of the electron-positron production rate using neutrinos in an 
intense background magnetic field. The computation is done for the process v — > vee (where v 
can be v e , z/ M , or v T ) within the framework of the Standard Model. Results are given for various 
combinations of Landau-levels over a range of possible incoming neutrino energies and magnetic 
field strengths. 
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I. INTRODUCTION 



Neutrino interactions are of great importance in astrophysics because of their capacity to serve as mediators for 
the transport and loss of energy. Their low mass and weak couplings make neutrinos ideal candidates for this role. 
Therefore, the rates of neutrino interactions are integral in the evolution of all stars, particularly the collapse and 
i-C ' subsequent explosion of supernovae, where the overwhelming majority of gravitational energy lost is radiated away in 
the form of neutrinos. 

Neutrinos have held a prominent place in models of stellar collapse ever since Gamow and Schoenberg suggested 
their role in 1941 While supernova models have progressed a great deal in the last 65 years, the precise mechanism 
i-C . for explosion is still uncertain. A common feature, however, among all models is the sensitivity to neutrino transport. 
Neutrino processes once thought to be negligible now become relevant, and this has inspired many authors to calculate 
rates for neutrino interactions beyond that of the fundamental "Urea" processes 

> 

I/"") ' p e — > n v e 



n — > p e v e . 

Recent examples include neutrino-electron scattering, neutrino-nucleus inelastic scattering, and electron-positron pair 
annihilation 0, 0]. Furthermore, the large magnetic field strengths associated with supernovae (10 12 -10 17 G) are 



likely to cause significant changes in the behavior of neutrino transport. 

While the the electromagnetic field does not couple to the Standard Model neutrino, it does affect neutrino physics 
by altering the behavior of any charged particles, real or virtual, with which the neutrino may interact. A number 
of authors have considered such effects on Urca-type processes @, H, @, 0| and on neutrino absorption by nucleons 
(and its reversed processes) H, Furthermore, Bhattacharya and Pal have prepared a very nice review of other 
j_J ■ processes involving neutrinos that are affected by the presence of a magnetic field [111 ]. 

The problem of interest in this work is the production of electron-positron pairs with neutrinos in an intense 
magnetic field 

v —>■ v e e. (1) 

Normally this process is kinematically forbidden, but the presence of the magnetic field changes the energy balance 
of the process, thereby permitting the interaction. 

Stimulation of this process with high-intensity laser fields has been shown to have an unacceptably low rate of 
production [l^ |. but such an interaction could have important consequences in astrophysical phenomena where large 
magnetic field strengths exist. The process would most likely serve to transfer energy in core-collapse supernovae [l3(. 
However, Gvozdev et al. have proposed that its role in magnetars could even help to explain observed gamma-ray 
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(a)Neutral current reaction 



(b)Charged current reaction 



FIG. 1: Possible diagrams considered for the process v — » vee. Both diagrams contribute for electron-type neutrinos, but 
only the neutral current reaction (FIG. 1(a) I contributes for and v T . 



bursts fl4j ]. The interest in this reaction has led to a previous treatment in the literature [Ta ]. but those authors 
present results for two special limiting cases: (1) when the generalized magnetic field strength eB is greater than the 
square of the initial neutrino energy E 2 , and (2) when the square of the initial neutrino energy E 2 is much greater than 
the generalized magnetic field strength eB. In both cases the incoming neutrino energy E is much greater than the 
electron's rest energy m e . In this paper we present a more complete calculation of the production rate as mediated by 
the neutral and the charged-current processes (FIG. [J). We present the results of the calculation for varying Landau 
levels, neutrino energies, and magnetic field strengths. A comparison with the approximate method is also discussed. 



II. FIELD OPERATOR SOLUTIONS 



As we have pointed out in section [H the standard model neutrino can only be affected by the electromagnetic field 
through its interactions with charged particles. This means that for the process v — ► vee the Dirac field solution for 
the final state electron and positron must change relative to their free-field solutions. The magnetic field will also 
change the form of the IF-boson's field solution which can mediate the process when electron neutrinos are considered. 
However, in our analysis we take the limit that the momentum transfer for this reaction is much less than the mass 
of the VF-boson (Q 2 <C tn^) and ignore any effects the magnetic field may have on this charged boson. Thus, in 
this section we review the results of our derivation of the Dirac field operator solutions for the electron and positron. 
We closely follow the conventions used by Bhattacharya and Pal and refer the reader to their work [9] for a detailed 
derivation. The reader who is familiar with these solutions may wish to begin with section IIIII where we calculate the 
production rate. 

We choose our magnetic field to lie along the positive z-axis 

B = B Q k (2) 
which allows us some freedom in the choice of vector potential A(x). We make the choice 

A*(a;) = (0,-yB,0,0) (3) 

both for its simplicity and its agreement with the choice found in reference Q . This choice in vector potential leads us 
to assume that all of the y space-time coordinate dependence is within the spinors. The absence of any y dependence 
in, for instance, the phase leads us to define a notation such that 

= (*,*,{),*) (4) 

and 

% = (V x ,0,V z ) , (5) 

where V is any 3- vector. 
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A. Electron field operator 



Solving the Dirac equation for our choice of vector potential (Eq. [3]) results in the following electron field operator 

d 2 p v lE n + m 



m*) -EE 



(2tt)2 



2E n 



u°(p v ,n, y) e-<™ a s e * n + v s (p^n, y) e +i ™ S^J. 



n=0 s=± • 

where the creation and annihilation operators obey the following anti-commutation relations 



(6) 



(7) 



In Eq. [5] we sum over all possible spins s and all Landau levels n where E n is the energy of fermion occupying the 
n th Landau level 



E n = y/p 2 z +ml + 2neB . 



n > 0. 



The Dirac bi-spinors are 



and 











(8) 



(9a) 



v + (Pv,n,y) 



V 



\/2neB j It 
E n +m e n ^~\ 





J 







(9b) 



The I m (0 are functions of the Hermite polynomials 



'eB 



1/2 



2™ m! x /tF 



where the dimensionless parameter £ is defined by 



e± = VeBy± 



'eB 



(10) 



(11) 



Recall that the Hermite polynomials H m (^) are only defined for nonnegative values of m. Therefore, we must define 
2— 1(£) = 0. This means that the electron in the lowest Landau energy level (n = 0) cannot exist in spin-up state and 
the positron in the lowest Landau energy level cannot exist in the spin-down state. 

The normalization in Eq. (fT0|) has been chosen such that the functions I m (0 obey the following delta- function 
representation [161 p. 86] 



S(y - y') 



\dy/dt\ 



'eB 



n=0 



5(y-y') = £/n(02»(O- 



(12) 



For convenience we choose to normalize our 1-particle states in a "box" with dimensions L x L y L z — V such that 
the states are defined as 



|e) = |pi^,ni,si) 



1 



(13a) 
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|e) = \p2 V ,n 2 ,s 2 ) = — L= S^„ 2 |0) , (13b) 



and the completeness relation for the states is 



n=0 s=± 



1 = Y.J2 I 4^^*IA»n»*><PV.n.*l ■ ( 14 ) 



B. Spin sums 

In order to evaluate the production rate for our process, we must derive the completeness relations for summations 
over the spin of the fermions. For a detailed calculation of the rules see reference [l(|. The results of the calculation 
are as follows 

J2 u s (p v ,n,y')u s (p v ,n,y) = (2(E n + m^y 1 ( [m e (l - a 3 )+ + ^ 7 5 ] /«(£_)£.(£-) 
«=+,- ^ 

+ [m e (l + <7 3 )+ ^|| 7 5 ] J„_i(£_)J»_i(£_) 

+ ^2neB (7 1 + ^7 2 ) /^(^ )!„(£_) 



+ V2^B( 1 1 -i^LaiOln-iit-)} (15a) 



and 



^^(^,n,yK(^,n,yO - (2(E n + m e )) l { [-m e (l - a 3 )+ ^+ ^ 7 5 ] /„(£+)/„(&) 

a=± 



+ [-m e (l + a 3 )+ ^||- ^|| 7 5 ] 7„_ 1 K+)7„_ 1 (^) 
- (7 1 + l7 2 ) 7„_ 1 (e+)J„(^) 

-V2^B( 7 1 - l 7 2 )/„(e+)/„-i(^)j, (15b) 



where 

Pl = (£,0,0, p.) (16) 
gjf = (p„0,0,£). (17) 

The above results have been derived using the standard "Bjorken and Drell" representation for the 7-matrices [l7j 



C. Neutrino field operator 



Having no charge, the neutrino's field operator solution ^(x) is not modified due to the magnetic field. Wc present 
it here for easy reference 

M*) = E / j^f [ uS w e ~ ip ' x a - + wS (p) e+ip ' x ' ( 19 ) 

where the creation and annihilation operators obey the conventional anticommutation relations 

{a*,aj't} = = (2nf 6 ss 'S 3 (p-f). (20) 
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The neutrino bi-spinors follow the standard spin sum rules 



(21) 



s=± 



s=± 



where we take the Standard Model neutrino mass to be zero. 

With "box" normalization the 1-particle states for the neutrino are 



satisfying the completeness relation 



1 



d 3 p 



V \p, s) (p, s\ 



(22) 



(23) 



III. THE PRODUCTION RATE 

The quantity of interest for the process v — > vee in a background magnetic field is the rate at which the electron- 
positron pairs are produced T. The production rate is defined as the probability per unit time for creation of pairs 



V 

r = lim — 



(24) 



where T is the timescale on which the process is normalized. We begin by finding the probability V of our reaction 



V 



^ J (2tt) 3 J (2tt) 

l,n2=0*' ^ ' J ^ ' 



d 2 p\ v 



2 L X L Z 



d 2 p 2# 
(2tt) 2 



L X L Z Y Up 1 > s '; Ply 



,n 1 ,Si 1 p 2v ,n 2 ,S2 



S 



p,s 



(25) 



S,S' ,Si ,S2 



In Eq. [2U quantities with the index 1 correspond to the electron, those with index 2 to the positron, the primed 
quantities to the final neutrino, and the unprimed quantities correspond to the initial neutrino. 



A. The scattering matrix 



The scattering matrix 



S=(p ,s ;pi v ,nx,si;p2y,ri2,S2 



S 



p, s 



(26) 



naturally depends on the flavor of the neutrino. While the process involving the electron neutrino can advance through 
either the charged (W) or neutral (Z) current, the muon (or tau) neutrino can only proceed through the latter. For 
this reason we will break the scattering matrix into a neutral component 



Sz = \P ,s ;pi v ,rii,si;p2y,ri2,S2 



S 2 



and a charged component 



Sw = [P ,s \pi v ,ni,si 1 p2 V ,n 2 ,s 2 



S, 



p,s 



p,s 



(27a) 



(27b) 



where the scattering operators are defined by the Standard Model Lagrangian as 



Sz 



2 3 cos 2 Ow sin 2 9w 



d^^rigv-gl^M^z^x) / dV^'h CT (l- 7 5 )^(z')^') (28a) 



'ir 



2 3 sin 2 6 W 



d 4 xMx)r (1 - 7 5 ) 1> v .[x)W-{x) / dVtJ/) 7 (1 - 7 5 ) M*)W+&) > 



(28b) 
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and 9yy is the weak-mixing angle, v\ indicates a neutrino of any flavor, v e refers to a electron neutrino, and the vector 
and axial vector couplings for the electron are 

9v = -i + 2sin 2 0w (29a) 
9a = ~\- (29b) 

In our analysis we will be using incoming neutrino energies that are well below the rest energies of the Z and W 
bosons. Therefore, we can safely make the 4-fermion effective coupling approximation to the Z and W propagators 

{Q\T{Z^x)Z a (x'))\0) -» 8\x-x')^ (30a) 

m z 

(0\T(W-(x)W+(x'))\0) - 5 4 (x-x')^f. (30b) 

m w 

After making this approximation our expressions for the scattering operators simplify to 

s z = % ! d^^rigv-gli^M^^hA 1 -^ 5 )^) ( 31a ) 



S w = % I d*xil> e (xW(l-^)rl> v .(x)il> v .(xftp(l-^)M*)> (31b) 



V2 
Gp 

where Gf/V2 — e 2 /(8 sin 2 6w m \y)> an d we have made use of the fact that cos 2 9w = mw/ m %- 

After substituting of the scattering operators (Eqs. (l3Tj) ) into the expressions for the components of the scattering 
matrix (Eqs. (f2"T)0 . we can use our results from sections III Al and III CI to write the components in the form of 



i(2n) 3 5 3 (p g -p' v - P v a ~P V ,2) 
L x L Z V 



&Z/W — j — J— 7} -M-z/W > \ 61 ) 



where 



-iG F (E ni + m e )(E n2 + m e ) _ / , 5 , 

MZ = M V EWK-^ 2 « (p)7„(1-7 6 )«(p) 

x j dye^y-P> u s i (p^, ni ,y) ^ (g c v - g Al 5 ) v s * (p 2in n 2 , y) (33a) 

iG t 



M 



w 



\ {E ni +m e )(E n2 +m e ) _ s , , , 5 , 
V EE>E ni E n2 « (P)7„(l-7)« (P2 V ,n 2 ,y) 



2 2 V2 

x [ &ye i{ - p »- p >u si (p lv ,n u y)j tl (l- 1 5 )u s (p). (33b) 



The reversal of sign on Eq. (|33b[) relative to Eq. ((33a|) is from the anticommutation of the field operators. The 
scattering amplitude for the charged component Mw can be transformed into the form of the neutral component 
M z by making use of a Fierz rearrangement formula 



"i7m (! - 7 5 ) U2 "37^ (1 - 7 5 ) u 4 = -ui7 M i 1 " 7 5 ) "4 u 3 ^ (l - 7 5 ) u 2 , (34) 



such that 



hA -iG F (E ni +m e )(E n2 +m e ) , . 5 , 

Mw = WfV — e^eZeZ 2 u (pMi-tVM 



x /dye^- p «) ! '^ 1 (p ll ,,n 1 ,y)7^(l- 7 5 )^(p 2l ,,n 2 , 2 ;) . (35) 



With the rearrangement of Mw m Eq. (|35p . we can now express the scattering amplitude in terms of the type of 
incoming neutrino. The muon neutrino can only proceed through exchange of a Z-boson, so its scattering amplitude 
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is just that of Mz 

M Va = M z 



iG F (E ni + m e )(E n2 + m e ) _ s > , , 5 , 

* = ^vlV ewKJ^ 2 u <pMi-tV(p) 

x / dye^y-P>u^ (pi v ,ni,y) 7 M (Gy - 7 s ) v S2 (p 2v ,n 2 ,y) . (36) 



The scattering matrix for a tau neutrino, and the subsequent decay rate, is exactly the same as the muon neutrino. 
We will keep the notation as for simplicity. 

The electron neutrino has both a Z-boson exchange component and an VF-boson exchange component. Therefore 
we must add the amplitudes to find its scattering amplitude 



M„ = M z +M 



-iG F (E ni + m e )(E n2 + m e ) , , , 5N 



x dye^y- p >u^(p lv ,n 1 ,y)^(G v - 7 5 )v S2 (p 2v ,n 2 ,y) . (37) 



Note that the scattering amplitudes for electron (Eq. l37|) and non-electron neutrinos (Eq. l36|) depend on a generalized 
vector coupling Gy defined by 

G v = 1 ± 4 sin 2 9 W . (38) 

We see that the scattering amplitudes for an incoming electron neutrino versus an incoming muon neutrino differ 
only in the value of the generalized vector coupling and an overall sign. And the overall sign will be rendered 
meaningless once the amplitude is squared. Therefore, we choose to make no distinction between the two processes, 
other than keeping the generalized vector coupling as Gy , until we discuss the results in section IIVI 



B. The form of the production rate 



Having determined the scattering matrix S and scattering amplitude M in section Till A[ we can now make series 
of substitutions of those results to find the expression for the production rate T. We begin by substituting the form 
of the scattering matrix (Eq. (f3"2")) ) into the expression for the production rate (Eq. 



lim — 



= lim T- 1 V 

T,V^oo ^— ' 
00 

= lim T- 1 V 



ni ,712—0 



dV 

(2tt) 3 

dV 

(2tt) 3 



V 



V 



(27T) 2 

d 2 pi # 

(2tt) 



L X L Z 



2 L X L Z 



(27T) 

d 2 p 2# 
(2tt)2 



;pi v ,n 1 ,si;p 2v ,n 2 ,s 2 



S 



p,s 



L X L Z 



i(27r) 3 J 3 -Py - p^i - ^,2) 



L x L z V 



M 



r = T jim co (2^T^)- 1 f; | dV | d 2 p lt , | d 2 ^ (<5 3 (^ -P' v - V V ,1 - Pv.2)) 2 \M\ 2 , 



m ,n2—0 



where |A4| is the square of the scattering amplitude after summing over spins 



\m\ 2 = y, \ M \ 2 - 



(39) 



(40) 



S,S' ,Si ,S2 



We can simplify the square of the 3-dimensional delta function by expressing one of the 3-dimensional delta functions 
as a series of integrals over space-time coordinates 



d 3 y 
(2tt) 3 



(41) 



By using the remaining set of delta functions to reduce the exponential to unity, we can write the integrand in 
of the dimensions of our normalization "box" 



P v , 2 )) 2 = 5 3 {p v -p v -p v j-p va ) 



(5 3 ( Plf -p' v - p v ,i - p v , 2 )) = 5 3 (p v -p' v - p v ,i - P v a) 



d 3 y 

(2tt) 3 
TL X L Z 
(2tt) 3 ' 



With the above result for the square of the delta function, the production rate in Eq. (|39|) simplifies to 



y rii,Ti 2 =0 

The square of the scattering amplitude goes as the product of two traces 

S,S f ,S1 ,S2 

_ (^f (E ni +m e )(E n2 +m e ) 
2 7 EE'E ni E„ 2 

x^5 s (p) 7ff (l- 7 5 ) u s \p')u s ' (j>') lf , (l- 7 5 )u s (p) J dye l< - p y- p > J dy' e' l{p y- p 'y )v ' 

s.s' 

x E v S2 (p2v, n 2 ,y') 7 CT (G± - 7 5 ) u Sl (p lv ,m,y') u Sl (f llt ,m,y) 7" (G± - 7 5 ) v S2 (p 2v ,n 2 , y) 



\M\ 2 = Of (EE'E ni E n2 ) 1 Jdye^y- p 'y^ J dy' e^-^Tr { 7ff (l - 7 5 ) / 7 m i 1 - 7 5 ) P'} 
xTr( 7 ff (G±- 7 5 ) \{m e (l-a 3 )+ ^||+ ^i||7 5 ) ^(C-.x^ni^-.i) 



+ V2^B (7 1 + i 7 2 ) I^xiC^Indt-.i) 
+ ^2n ie B (7 1 - ij 2 ) I ni (C,l)Ini-l{Z-,l) 
+ (m e (l + a 3 )+ ^l|7 5 ) 

x r (G± - 7 5 ) f(-m e (l - <7 3 )+ A||+ ^ 2 || 7 5 ) /„ 2 (^, 2 )/„ 2 (e+,2) 



- V2^B (7 1 + Z7 2 ) /„ a -i(^, 2 )/„ 2 (e+,2) 

- ^2^8 (7 1 - ^ 7 2 ) /„ 2 (^, 2 )/„ 2 -i(e+, 2 ) 

+ (-m e (l+<7 3 )+ ^ 2 ||7 5 ) J„ 3 -i(eV,a) J n a -i(f+ J a) 



where we have used our result for the summations over spin from Eqs. (|15p and (12 1[) . 
The space-time dependence of Eq. (|44l) can be factored into terms like 



dye^- p > /„(£_,!) 7 m (£+, 2 ) 



and 



dj/e-^-*'i)*'j n (^_ il )J ro (r + , 3 ), 



where the /„ jm are functions of the momenta in the problem. 
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We have included a detailed calculation for the general form of I n . m in appendix [XI but we only present the result 
here 

I = ' /ie"" 2 /2 e #o {Vx + ir)y) m-n L rn-n ^2) m > n >Q 

M^ v2/2 e# ° i-Vx + iVv) n ~ m L?n m (V 2 ) , n > m > o' 



where 



Plx +P2x 



Py^Py 



(48) 



% = 1m (49) 



(Py-Py)(Pi -Vi) 
2eB 



(50) 



V 2 = Vl+V 2 y , (51) 

and L™~ n {r} 2 ) are the associated Laguerre polynomials. 

The full results of the traces and their subsequent contraction are nontrivial but have been included in appendix[EJ It 
is important to note, however, that the only dependence on the ^-components of the electron and positron momentum 
is that which appears in Eq. (|47p for I n ^ m . Furthermore, we notice that all terms in the averaged square of the scattering 
amplitude have factors that go as a product of I nm and /*, m ,. Therefore, the coefficient e i( ™ in Eq. (|47[) will vanish 
when this product is taken. The only remaining x-dependence of these two momenta appear as their sum in the 
parameter r\ x = (jp\ x + P2x)/V2eB. This helps to simplify the phase-space integral for our production rate (Eq. (|43P) 
which is proportional to 

Tex lim — / dpi x [dp 2 x- (52) 

L y >00 Ly J J 

If we make a change of variable from the a;-component of the positron momentum p 2x to the parameter r) x , the 
relationship in Eq. (|52p is rewritten as 

T oc lim ^ 2eB f dp lx [ drj x . (53) 

L y ^co L y J J 

Because there is no longer any explicit dependence on the ^-component of the electron's momentum p\^ x in the 
averaged square of our scattering amplitude, we can simply evaluate the integral 

To evaluate this integral we must determine its limits. As discussed previously, we have elected to use "box" normal- 
ization on our states. This means that our particle is confined to a large box with dimensions L x , L y , and L z . The 
careful reader will note that we have already taken the limit that these dimensions go 

to infinity in some places, particularly in Eq. (|A5[) . but it is imperative that we be cautious here, as we could naively 
evaluate the integral over p\^ x to be infinite. 

Physically, the charged particles in our final state act as harmonic oscillators circling about the magnetic field lines. 
While they are free to slide about the lines along the z-axis, the particles are confined to circular orbits in the x and 
y-directions no larger than the dimensions of the box. For a charged particle undergoing circular motion in a constant 
magnetic field, the x-component of momentum is related to the y-position vector by 

p x = ~eQBy (54) 

where Q is the charge of the particle in units of the proton charge e — |e|. Therefore, the limits onpi^ are proportional 
to the limits on the size of our box in the y-direction. The integral over the electron's momentum in the x-direction is 

r eBL y /2 

dpx j x = eBLy , (55) 

-eBLy/2 

and the result helps to cancel the factor of L y that already appears in the form of the production rate. We can now 
safely take the limit that our box has infinite size, and the production rate now has the form 



r = £ f dV f dp lz I dp 2z f d Vx V2^B S 3 ( Pv -p' v - p vA - ft, i2 ) eB ^} . (56) 
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IV. RESULTS 



In our expression for the total production rate (Eq. (|56|) ), one will notice is that there is a sum over all possible 
values of the Landau levels. As a consequence of energy conservation, upper limits do exist for the summation over 
the electron's Landau level n\ 



E — E + E ni + E ri2 
E > E n , + m P 



E — m e > \Jm\ + 2nieB 
5-2 
2eB 



E(E - 2m e ) , x 

ni < 1 „„ D , (57) 



and a similar one for the positron's Landau level 

E = E + E ni + E„ 



E > y/m 2 e + 2meB + E n 



E - y/ml + 2nieB > yjm 2 + 2n 2 eB 



n, < ± ^-^ . (58) 

These relationships help to constrain the extent of the summations. Physically, these constraints can be thought of 
as limits on the size of the electron's (or positron's) effective mass, where the electron (or positron) occupying the n th 
Landau level has an effective mass 



= y rn 2 + 2neB (59) 

and energy 



En = VpI+™* 2 - (60) 

For low incoming neutrino energies and large magnetic field strengths (eB > m 2 ), the constraints put very tight 
bounds on the limits of the summations. However, higher incoming energies and low magnetic field strengths impose 
limits that still require a great deal of computation time. For instance, at threshold (E = 2m e ) there can exist only 
one possible configuration of Landau levels {n\ = n 2 = 0), while at an energy ten times that of threshold and a 
magnetic field equal to the critical field (B = B c = m 2 /e = 4.414 x 10 13 G) there are nearly 7000 possible states. At 
the same magnetic field but an energy that is 100 times that of threshold, there are almost 70 million states. However, 
for incoming neutrino energies less than a certain value 

E < m e + yjml + 2eB (61) 

only the lowest Landau level is occupied, ri\ , n,2 = 0. And even at energies above, yet near, this value we expect that 
production of electrons and positrons in the n-y , n 2 — level is still the dominant mode of production because it has 
more phase space available. 

Production rates at the 0, Landau level are presented in FIG. [2] for both the electron and muon neutrinos. (All 
of the results for muon-type neutrinos are valid for tau-type neutrinos.) One interesting feature of these results is 
the flattening out of the rates at higher energies. The energy region at which this flattening begins increases with 
increasing magnetic field strength, and it appears to be in the neighborhood of energies just above the limit set in 
Eq. (|6ip. At energies in this regime we expect that modes of production into other Landau levels are stimulated, 
which helps to explain why the behavior of the 0, production rates change above this area. 

We should note that the results given in this work are all for an incoming neutrino traveling transversely to the 
magnetic field. The rates are maximized in this case as can be seen in the example found in FIG. [3] for an initial 
electron neutrino with energy E Va — 20m e in a magnetic field equal to the critical field B = B c = m 2 /e. 

For comparison purposes, the production rates for other combinations of Landau levels have been calculated. These 
include the 1,0 and 0, 1 cases (FIG. gj}, the 20,0 and 0,20 cases (FIG. |5J, and the 10, 10 case (FIG. |BJ. The first 
noteworthy feature of these results is that the production rates are decreasing at higher Landau levels. Because the 
energy required to create the pair goes as 

-^pair — E ni -\- E n2 



lp\ z + 2n,ieS + m2 + ^J p 2 z + 2n 2 eB 
E pair > y/2 ni eB + m 2 + ^2n 2 eB + m 2 , 
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FIG. 2: Production rates for the m , 112 = Landau levels where T is the rate of production, E v is the energy of the incoming 
neutrino, and the magnetic field is measured relative to the critical field B c = 4.414 x 10 13 G. All plots are for a neutrino that 
is perpendicularly incident to the magnetic field. 



the available phase space for the process should decrease in the order 0,0 — > 0, 1 — * 0, 20 — * 10, 10. And as can 
be seen in FIGS. [_ QJ [_ and[_ the production rates fall off accordingly. 

Another interesting feature of these results is the apparent preference for the creation of electrons in the highest of 
the two Landau levels. That is, the rate of production is larger for the state ri\ = i , = than for n% = , = i 
(FIGS. H] and [5}. This behavior is especially significant over the range of incoming neutrino energies near its threshold 
value for creating pairs in the given states. Though the i,0 production rate is larger and increases more quickly 
in this "near-threshold" range than its 0, i counterpart, both curves plateau at higher energies, and their difference 
approaches zero. This difference is presumably caused by the positron having to share the Ws energy with the final 
electron-type neutrino. This also explains why such an effect is not seen for muon and tau-type neutrinos that only 
proceed through the neutral current reaction. 

It was mentioned in section U that previous authors have considered this process under two limiting cases [la ]. One 




FIG. 3: The production rate's dependance on the direction of the incoming neutrino.. The production rate is for the 0, Landau 
level with an electron of energy E v — 20m e traveling at an angle 9 relative to a magnetic field of strength equal to the critical 
field B = B c . Data is included for both an incoming electron- type neutrino (solid line) and a muon-type neutrino (dashed line). 
If we average over 8, then the average production rate is 1.38 x 10 -16 cm -1 for electron-type neutrinos or 2.94 x 10~ 17 cm" 1 
for muon or tau-type. 



is when the square of the energy of the initial-state neutrino and the magnetic field strength satisfy the conditions 
E v 3> eB rnr e . Under these conditions many possible Landau levels could be stimulated, offering a multitude of 
production modes. Therefore, it would be inappropriate to compare their expression to our results for a specific set 
of Landau levels. However, the second limiting case is for eB > E 2 m 2 .. This condition is slightly more restrictive 
than our condition for the energies below which only the lowest energy Landau levels are occupied (Eq. (|61|) ). In this 
regime our results for the 0,0 state are the total production rates, and we can compare our results to the expression 
derived by the previous authors [15| 

G 2 F (G± 2 + l) 

T = ^-5 J -eBE* (l + 0(E 2 /eB)) , (62) 

where we have taken the direction of the incoming neutrino to be perpendicular to the magnetic field's direction. 
Results of this comparison are shown in FIG. [JJ 

The results in FIG. [7] demonstrate the drawbacks of using the approximation in Eq. ((62|) . While the expression is 
very simple, it gives only reasonable agreement with the production rate at a magnetic field equal to 100 times that 
of the critical field (B = 100 B c ). Here it overestimates, at the very least, by a factor of two, and the inclusion of 
higher order corrections makes no significant improvement. One reason for the disagreement at this field strength is 
that there is only a very small range of energies that satisfy the condition eB > E? ^S> m 2 . Therefore at higher field 
strengths we should get better agreement, and we do. Closer inspection of FIG. [7] reveals that the differences are 
less than a factor of three for neutrino energies in the range 2 MeV < E v < 20 MeV, and the expression successfully 
provides a good order of magnitude estimation. Though the estimate will improve at higher magnetic field strengths, 
it begins to loose relevance as there are only a handful of known objects (namely magnetars) that can conceivably 
possess fields as high as 10 15 G. Even for these objects, fields stronger than 10 15 G cause instability in the star and 
the field begins to diminish [1 31 ] - 

Probing the limiting case E v \fe~B is imperative because our present work has already demonstrated nontrivial 
deviation from approximate methods for realistic astrophysical magnetic field strengths and neutrino energies near 
and below the value \J eB. But, as was mentioned previously, the number of Landau level states which contribute to 
the total production rate grows very rapidly in this higher energy regime, and we need to sum over these states. Future 
work will attempt to do these sums by using an approximation routine that can interpolate between rates for known 
sets of Landau levels. This will provide a flexible way to balance accuracy with computation time while determining 
when the production rate deviates from its limiting behavior. The significance of these deviations will only be known 
when a more complete understanding of the role that neutrino processes play in events such as supernova core-collapse 
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and in the formation of the resulting neutron star. This work aims to improve that understanding. 
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APPENDIX A: CALCULATION OF I n , m 

In section fill B I we discuss the fact that the squared scattering amplitude has coefficients that are integrals over the 
space-time coordinate y 

In,m = J dye^y~P>I n (^. 1 )I m ^ + . 2 ). (Al) 
In this appendix we will derive the result after integrating over y. 
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By defining new parameters 

C = V^By (A2) 

Q = p ix /V^B (A3) 

Co = (p y -p' v )/V^B (A4) 

and using the definition of £ (Eq. (jlip ) we can make a change of variable from y to £ and rewrite l n ,m as 



s lCoC /«(C-Ci)/ m (C + C 2 ) 



(A5) 



where the limits of integration are ±00 because we have taken the limit of L y as it approaches 00. The in 
Eq. (|10l) depend on the Hermite polynomials if„(£), which can be represented as a contour integral in the following 
way [3 Eq. (13.8)] 



-n—l —r+2tx 



(A6) 



Substituting this definition of the Hcrmitc polynomial into Eq. (fT0|) allows us to write the I n?m as 



dC ,, 



CoC 



'eB 



1/2 



'eB 



2 n n\ x/tt 



e-(^) 2 / 2 i/„(C-Ci) 



'eB 



1/2 



2 m to! \/tt 



e -(c+c 2 ) 2 / 2ifm(c + C2) 



(2"+ m n! to! ^)" 1/2 ^ dC e* Co c e -(f-^) 2 /2 H ^ C _ Cl ) e -(C+C 2 ) 2 /2 # m(c + ^ 



= (2 r ' 



,-1/2 



dCe 



<oC P -(C-Ci) J /2 p -(C+C2) 2 /2 



x /d^-n-ig-^+atCC-d)^ ^ dss - m -i e ^ +2s (c+c 2 ) 
2iri J 2ni 



(A7) 



Next, we isolate all of the £ dependence, interchange the order of the integrals, and perform the integration over £ 



(A8) 



Inti = / dC exp (-C 2 + C(Ci - (2 + <o + 2t + 2s)) 

J — OC 

= exp ((Ci - (2 + Ko + 2t + 2s) 2 /4) . 
Substitution of this result back into Eq. (|A7j) gives 

In,m = (2 n+m n\m^ 1/2 e -(K 1 +&) 2 +Co 2 )/4 e iCo(C 1 -6)/2 

x — I dtt-"- 1 e *(-Ci-6-HCo) ^ (/■dss-^e^+^+^e 23 * 

27TI / 27TI 



(A9) 



If to > n, then we can perform the integration over s first 

ml 
2ttI 



Int 2 = __ dss^V^+^+^V 3 * 



d" 



s(Cl+C2+iCo+2t) 



ds m 

(Ci + c 2 + *Co + 2ty 



s=0 



(A10) 



such that 



I n ,m = (2" +m n! m!) 1/2 exp (- ((Ci + C2) 2 + Co) /4) exp (i$,(£i - 0»)/2) 
n! / (Ci + C 2 + <o + 2t) m 

x 2^f ^ exp ( ( ~ Cl " (2 <o)) ■ 



(AH) 
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The integration over t is made easier by making the following changes of variable 

* - HP - w (A12) 

, Co (Ci - C2) (Py ~ Py)(Pi ~ P2) . . 

^ 2 - 2^B (AU) 

rf 1 = r] x ± ir^ (A15) 

7? 2 = 77+77- = 7/2+7/2 (A16) 

i = (u-77 + ) /V2. (A17) 
The integration over the variable t can now be written as 

l nt „ = n - I du (%/2u) m c ((u^+)/V2)(-V2 V -) 



2ni J (u- ?7+) n+1 

_ 2(™+ m )/ 2 g»7 _ u m e ~ u V~ 

du n 

= ri!2 ( ™ +m)/2 (?7 + ) (m -™ ) i™-"(77 2 ) 
where we have used the Rodrigues' representation for Laguerre polynomials pj, Eq. (13.47)] 

L k (x) = tJL— T— * n+k e~ x , n,k>0. (A18) 
71 w n\ dx n ' ~ v 



With the result from Eq. (lATSl) . we can now express the I n ,m as 

In,m = \/4 e_J72/2 fa* + 1 ~"V) - "1 > 71 > , (A19) 



m! 

For the case when n > m we first integrate over t in Eq. (|A9[) and follow a similar procedure to find 

' m! 



I 2 (-77, + i 7;,,)"-" 1 L n m m {if) , n > m > . (A20) 



APPENDIX B: RESULT OF TRACE 

We can express the trace result for the average of the squared scattering amplitude from from Eq. ((44)) as a sum of 
terms 

-,2 16 



= 2»EE G E ni E n S> T " < B1 > 



where the coefficients A4 depend on the products of I nm and I*, , defined in Eq. (147[) and presented in appendix [Al 
and the Tj are the parts that depend on the contraction of the traces in Eq. (I44|) . The results are as follows: 

Ai = Ini,n 2 Ini,n2 ( B2 ) 
Ti = Tr{ 7 -(G±- 7 5 ) [m(l - a 3 )+ + 4 nl 5 ] 7 M (G$ - 7 5 ) [-m(l - a 3 )+ ^ + ^| 7 5 ] } 
Ti = -2 7 (G± 2 - l) m 2 (p x p' x +p y p' y ) + 2 6 (G± + l) 2 (E - p z )(E' - p' z )(E ni + Plz )(E n2 + p 2z ) 

+2 6 (G±-l) 2 (E+p z )(E'+p' z )(E ni -pi z )(E n2 -p 2z ) (B3) 



^-2 — J -ni,n 2 -l 1 ni,n2 

T 2 = Tr{ 7 CT (G±- 7 5 ) [m(l-a 3 )+ ^i||7 5 ]7 M (G±- 7 5 )(-V2^B)(7 1 +H 2 )} 
T 2 = -2 6 (G± + l) 2 V2n 2 eB (p^ + i Py )(E' - p' z )(E ni + Plz ) 
-2 6 (G± - l) 2 V2^B (p' x + i + p z )(£? ni - Pl z ) 

t 3 = t; 

- In 1 ,n 2 -ll ni ,n 2 -l 

= Tr{ 7 CT (G±- 7 5 ) H1- ( t 3 )+^ 1 ||+^ 1 || 7 5 ] 7 ^(G±- 7 5 ) [-m(l + a 3 )+^ 2 ||- ^| 7 5 ]} 
= 2 7 (G± 2 - l) m 2 e (E+p z )(E'-p' z ) + 2 6 (G± + l) 2 (£7 + p z )(E' - p' z )(E ni +p lz )(E n2 - p 2 
+2 6 (G± - l) 2 (£+p z )(£;'-p'J(£ ni -pi z )(E n2 + P2z ) 

Ar — T T* 

^-5 — 1 n u n 2 1 n 1 ~l,n 2 

T 5 = Tr{ 7 ff (G± - 7 5 ) (T^B) + * 7 2 ) r (Gy - 7 5 ) [-m(l - <7 3 )+ jfe„+ ^ 2 ||7 5 ] } 
T 5 = 2 6 (G± + 1) 2 ^/2^B(p' x + ip' y )(E-p z )(E n2 +p 2z ) 

+2 6 (G± - l) 2 y/2meB ( Px + ip y )(E' +p' z )(E n2 -p 2z ) 

Aq = In 1 ,n 2 -ll ni -l. n2 

T 6 = Tr{ 7 CT (G±- 7 5 )(V2^)( 7 1 +H 2 )7 P (G^-7 5 )(-V / 2^)(7 1 +*7 2 )} 
T 6 = -2 6 (G v + l) 2 v/2nieB v^eB (p x +ip y ){p' x + ip' y ) 
-2 6 (G± - l) 2 v/^B VS^eB (p* + ip„)(p^ + ij/J 

^7 = ^ni,n 2 -d-l,n 2 -l 

T 7 = Tr{ 7 CT (G±- 7 5 )( v ^^)(7 1 +H 2 )7^G±- 7 5 )(- v ^B)( 7 1 -» 7 2 )} 
T 7 = -2 6 (G± + l) 2 y/2n x eB ^2n 2 eB ( Px - ip y ){p' x + ip' y ) 
-2 6 (G± - l) 2 ^2n 2 eB ( Px +ip y ){p' x - ip' y ) 

A J J* 

^8 L n\,n 2 — l J m- l,n 2 — 1 

T 8 = Tr{ 7 ff (G± - 7 5 ) (V2^B) (7 1 + ij 2 ) ^ (G v - I 5 ) [-m(l + ct 3 )+ jfe,,- ^| 7 5 ] } 
T 8 - 2 6 (G± + 1) 2 ^/2^B(p' x + ip' y )(E+p z )(E n2 -p 2z ) 

+2 6 (G± - l) 2 ^j2n x eB ( Px + i Py )(E' + P ' z )(E n2 +p 2z ) 
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a 9 = i ni -!, n2 r nun2 = a; (bis) 

T 9 = T* (B19) 

Aw = ^m-l,n2-l-d,n 2 = ^7 (B20) 

Tio = T 7 * (B21) 

An = ^m.na^ni.na-l = ^6 ( B22 ) 

Tn = T 6 * (B23) 

^12 = ^ni-l,n 2 -l-C 1)n2 _l = ^8 (B24) 

T12 = T* (B25) 

A 13 = ^ni-l.na-Ci-l.na ( B26 ) 

T 3 - Tr{ 7 CT (G±- 7 5 ) [m(l + a 3 )+ Ay- ^|| 7 5 ] ^ (G± - 7 5 ) [-m(l - a 3 ) + ^ + ^| 7 5 ]} 

T 13 = 2 7 (G± 2 - l) m^S - + p' z ) + 2 6 (G± + l) 2 (E - p z )(E' + p' z )(E ni - Pl z )(E n2 + p 2z ) 

+2 e (G$-l) 2 (E- Pz )(E'+ P ' z )(E ni + Plz )(E n2 - P2z ) (B27) 

A\i = ^m-l,n2-l^ni-l,n 2 (B28) 

r w - Tr{ 7 ff (G±- 7 5 ) [ TO (l + CT 3) +J J 1|r ^ 1||7 6]y'(G^-7 6 )(- v ^B)(7 1 +i7 2 )} 
T 14 - -2 6 (G± + l) 2 {^/2~^B){p x + i Pv )(E' +p' z )(E ni - Plz ) 

-2 6 (G± - l) 2 {y/2n 2 eB)(p' x + ip' y )(E - Pz )(E ni + Plz ) (B29) 

-4-15 = Jni-l,n2-fni-l,n 2 -l = ^14 ( B30 ) 

Tib = T; 4 (B31) 

-4l6 = ^ni-l,n2-l-d-l,n 2 -l (B32) 

Tie = Tr{ 7 CT (G±- 7 5 ) [m(l + ct 3 )+ ^y- ^|| 7 5 ] ^ (G± - 7 5 ) [-m(l + a 3 )+ jfe,,- ^| 7 5 ]} 
Tie - -2 7 (G± 2 - l) m^^+PX) + 2 " ( G v + (E + p z )(E' + p' z )(E ni - Plz )(E n2 -p 2z ) 

+2 6 (G± - l) 2 (E - p z )(E' - p' z )(E ni + Plz )(E n2 + P2z ). (B33) 



